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Abstract 

We will study the least square estimator 9t,s for the drift parameter 9 of the frac- 
tional Ornstein-Uhlenbeck sheet which is defined as the solution of the Langevin equa- 
tion 



Xf.s — 



off Xy,udvdu + B"f, {t,s) £[0,T]x[0,S] 
Jo Jo 



driven by the fractional Brownian sheet B°''^ with Hurst parameters a,f3 in (^,|)- 
Using the properties of multiple Wiener-Ito integrals we prove that the estimator is 
strongly consistent for the parameter 6. In contrast to the one-dimensional case, the 
estimator 9t,s is not asymptotically normal. 

2010 AMS Classification Numbers: 60G15, G0H07, 60G35, 62M40. 

Key words: fractional Brownian sheet, parameter estimation, multiple Wiener-Ito 
integrals, strong consistency. 



1 Introduction 

We consider the two-parameter fractional Ornstein-Uhlenbeck process defined as the solu- 
tion of the stochastic equation 
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Xt,s = -e I I X.^udvdu + B'^f, (t, s) E [0, T] x [0, 5]. (1) 
Jo Jo 

Here B'^'^ denotes a fractional Brownian sheet with Hurst parameters a, (3 E !)• We also 
suppose that Xq^ = X^ q = Xq ^ = for every t, s. Our goal is to estimate the unknown 
parameter 6 from the continuous time observation of the solution (^t,s){t,s)G[o,T]x[o,S]- 

The development of the stochastic analysis for fractional Brownian motion (fBm) 
naturally led to the study of the statistical inference for stochastic equations driven by this 
process. There already exists an important literature related to these aspects. We refer, 
among others to [3], [1], [9], [13], [22]. Statistical analysis of the stochastic differential 
equations (SDE) driven by the fractional Brownian sheet has been less considered. We refer 
to the paper |18] for the study of the maximum likelihood estimator for a SDE with additive 
fractional Brownian sheet noise (see also [6] or [T] for the case when the noise is a standard 
Brownian sheet). 

In this paper we propose a least square estimator for the unknown parameter 6 
following the approach in [9]. This estimator is obtained by formally minimizing with 
respect to 6 the expression 



T r-S 

Jo 



Xt^S + GXt^; 



dtds 



2 

dsdt. 



We obtain the following estimator 



io lo Xt,sdXt,s 



7T,S 



loIoXlsdtds 



(2) 



The integral with respect to dXt^s is understood as the sum of the standard Lebesgue integral 
—6 Jq Jq dtdsX^g and of the stochastic integral Xt^sdB^f which is a divergence type 

integral with respect to the fractional Brownian sheet S"'^ (it will be defined in Section 2, 
we also refer to [10], [TT], [I2j, [20], [21] for the stochastic integration with respect to B°''^). 
Using ([1]) and ([2]) we can write 



lo lo ^ty'l'^^t,s 

lolo^ldtds 



o a _ JO JO ""i,^— t,s 
Ot,S -0 T r-S y2 ■ • 



We will study the asymptotic behavior of the least square estimator Ot,s as T, S* — oo. 
Our tools are the multiple stochastic integrals and the Malliavin calculus. Actually, the 
nominator and the denominator of the right hand side of ([3]) can be expressed as multiple 
integrals of order 2 with respect to the fractional Brownian sheet and from this we will 
obtain concrete estimates for their moments. We will prove that the estimator ([2]) is a 
strongly consistent estimator in the sense that it converges almost surely to the true value 
of the parameter 9. This result is similar to the one-dimensional case (see [9|), however the 
approach presented in [9] is not possible to be followed for the two-parameter case, instead 
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we use among other tools, the hypercontractivity of multiples integrals. By contrary, in the 
two-parameter case, the least square estimator does not preserve the asymptotic normality 
as in the one-parameter case. This will be noticed at the end of our work by using criteria for 
the asymptotic normality of sequences of multiple integrals in terms of Malliavin calculus. 

Our paper is structured as follows. Section 2 contains some preliminaries on multiple 
integrals and fractional Brownian sheet. In Section 3 we discuss the relation between the 
solution to ([1]) and the Bessel function of order 0. Section 4 contains the proof of the 
consistency of the least square estimator while Section 5 is devoted to a discussion about 
the asymptotic normality of the estimator. 

2 Preliminaries 

Let us introduce the elements from stochastic analysis that we will need in the paper. 
Consider % a real separable Hilbert space and {B{ip), ip € H) an isonormal Gaussian process 
on a probability space (fi, A, P), that is, a centered Gaussian family of random variables such 
that E {B{ip)B{tp)) = {ip,ip)'}{. Denote by In the multiple stochastic integral with respect 
to B (see [16]). This In is actually an isometry between the Hilbert space 7^®" (symmetric 
tensor product) equipped with the scaled norm ^/riJ.\\ ■ ||-^8)n and the Wiener chaos of order 
n which is defined as the closed linear span of the random variables Hn{B(ip)) where <p S 
H, W^Wu = 1 and Hn is the Hermite polynomial of degree n > 1 

«„M = tirexp(f)|;,(exp(4)), ..M. 

The isometry of multiple integrals can be written as: for m, n positive integers, 
^{In{f)Ini{g)) = n\{f,g)-y^s,n if m = n, 

E(/n(/)/™(<7)) = if m^n. (4) 

It also holds that 

W) = In{f) 

where / denotes the symmetrization of / defined by 

f{xi, ...,Xx) = ^^ ^ fiXail),- ■ ■,X„(^n))- 

The Malliavin derivative acts on multiple integrals F = In{f) in the following way: for 
every s 

DJn = n/„_i(/(-,s)) 

where " •" above denotes n — 1 variables. We recall the following hypercontractivity property 
for the L'P norm of a multiple stochastic integral (see p4l Theorem 4.1]) 

E|I^(/)|'"<c™(EI^(/)2)" (5) 
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where Cm is an explicit positive constant and / G T-L®"^. 

In this work we use Malliavin calculus and multiple integrals with respect to the 
fractional Brownian sheet (fBs). Let us define this process and its associated Hilbert space. 
The fBs with Hurst parameters a,/3 E (0, 1), (Sf/, s G [0, T] x [0,5]) is a zero mean 
Gaussian process with covariance 

given for all t,u£ [0,T]2 and s,v £ [0,S]^. 

We assume that B'^'f^ is defined on a complete probability space (O,^, P) such that 
A is generated by B°''^. Fix a time interval [0, T] x [0, S], denote by ^ the set of real valued 
step functions on [0, T] x [0, S] and let n°''(^ be the Hilbert space defined as the closure of 
with respect to the scalar product 

(l[o,t]x[o,s]) l[o,M]x[o,t)])w°"3 = T^"'{t,u)TZ'^{s,v) 
where 7^"(^, u)7^^(s ,v) is the covariance function of the fBs, given in 1^. The mapping 
l[o,t]x[o,s] ' — ^ ^t, f be extended to a linear isometry between 'H°''^ and the Gaussian 
space T-L^'^ spanned by B"''^ which is a closed subspace of L'^{Q, A,F). We denote this 
isometry by (p i — y B°''^{(p). Fix a, /3 > ^, in this case we have that for every f,g £ H"^'^ 
the scalar product has the form 

S T' S 

{f,g)^c,,0 = c{a)c{p) I I I [ fia,b)gim,n)\a-m\'^''-^\b-n\'^^-^dadbdmdn (7) 
Jo Jo Jo Jo 

and c(a) = a(2a — 1). 

3 About the solution 

The equation (JT]) has been studied in several papers (see [8], [E]). It has been showed 
that for 9 > and a, /3 > ^ equation ([1]) admits an unique strong solution which can be 
expressed as 

Xt,s= I I f{t,s,to,so)dB^^i^ (8) 
Jo Jo 

where 

f{t,s,to,so) = l[o,f](to)l[o,.](^o) ~ ~ (9) 

n>0 ^^'^ 

We will call the process X solution to ([T]) as the fractional Ornstein-Uhlenbeck sheet. It 
is a Gaussian process since it is given by a multiple integral of order 1 (Wiener integral 
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actually) with respect to the Gaussian process S"'^. We mention that the solution to 
([1]) behaves differently as its one-dimensional counterpart which is fractional Ornstein- 
Uhlenbeck process introduced in [5], this will make our analysis quite different from p]. For 
example, we note that the solution of some stochastic differential equations driven by the 
Brownian sheet or fractional Brownian sheet, which are positive in the one-parameter case, 
can take negative values with strictly positive probability (see [15] or [16j). 

A key element of our analysis is the fact that the solution X (more precisely the 
kernel / of the solution) can be expressed is terms of the Bessel function of the first kind. 
Let us consider the Bessel function of order given, for every x E M, by 

n>0 

This Bessel function admits the integral representation, for every x G M 

1 r 

Jq{x) = — / cos {xsinp) dp. 
Jo 

The kernel / in Q of the solution (^t,s)t,se[o,T]x[o,5] can be expressed as 

f{t,s,u,v) = l[o,t](n)l[o,s](^)Jo (2yje{t - u){s - v)) (10) 

= Mo,t]{u)l[o,s]iv)'^ COS (2y/9{t - u){s - v) sinp^ dp. 

Let us also recall the following property of the Bessel function (see e.g. p]) which will play 
an important role for our estimates: for x large enough 

/ 2 vr 

Mx) \ — cos(2;--) (11) 
V vrx 4 

(the symbol ~ means that the two sides have the same limit as x — )• (X)). 

4 Asymptotic behavior of the least square estimator 

In this section we study the asymptotic behavior of the estimator 6t,s defined in ([2]). More 
precisely, we will show that this estimator is strongly consistent for the parameter 9, that 
is, Ot,s converges to 6 almost surely as T, S" — t- oo. To this end we will analyze separately 
the nominator and the denominator appearing in the right hand side of the expression ([3]). 
Let us start with the study of the nominator. It can be written as the double stochastic 
integral 



/ / Xt,sdB^f ■.= FT,s:=h{f{u,v,t,s)) 
Jo Jo 



(12) 



where the kernel / is given by ([9]) and the integral I2 acts with respect to the variables 
{u,v), (t,s). 

We will estimate first the norm of Fts- We have the following result. 
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Proposition 1 For every e > and for a,f3^ |), 

E ^^-2a+^-eg-20+l-e Xt^sdB^f'^ ^ when T,S^oo. (13) 

Moreover for T, S large enough we have 

rT rS 

lo Jo 

where C is a strictly positive constant not depending on T, S. 

Proof: We calculate the norm of the random variable I2 {fiu, v, t, s)). By the isometry 
property of multiple integrals (jl]) and since ||/||c^c«,/3)(gi2 < ||/||(-^a,/3)®2 this norm can be 
handles as follows 



E ( r-2"+3 5-2^+3 / I Xt sdBff ] < C 



It.s < / dtdtoduduQ / dsdsodvdvQ 

x/(i, s, u, v)f{to, so, uo, vo)\u - uol'"-'!?; - ?;o|'^"'|t - to|'°"'|s - sol'^"' 

T I't i-T i-to rS t-s i-S t-SQ 

dt du dto / duo / ds dv dsQ / dvQ 
Jo Jo Jo Jo Jo Jo Jo 

X Jo (2y/e{t-u){s-v)'^ Jo (2V^(to-no)(so-^^o)) 
x|t - toP"-> - no|2"-2|5 - 5oP^-> - vo\^^-\ 

By making the change of variables i = u = ^ and similarly for the other variables, we 
obtain 

/■I ft rl rto rl rs rl /-sq 

Its = T^"-4r4^4/3-4^4 / dto / duo ds dv dso / dvo 

Jo Jo Jo Jo Jo Jo Jo Jo 

X Jo (2^9{t-u){s-v)Ts') Jo (2^/e{to-uo){so-vo)Ts) 

x\t- toP"-'|n - 7XoP"-'|s - sof'-^lv - vo\^^-' 
._ t^»s^^Ut,s- (14) 



Using the asymptotic behavior of the Bessel function (jllh . we have that 

Jo (2^ye{t-u){s -v)TS^ 



1 , -^T,S~ 

{TS)-4+^ 



for almost every t,u,s,v £ (0,1) and for every e > 0. We will next apply the dominated 
convergence theorem. To this end, using again relation (|11|) . it suffices to show that the 
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integral 



1 i-t i-l i-Ui j-l r-s r-l /■So 

dt I du I dto / duo / ds dv dsQ / dvo 



Jo JO JO JO JO JO JO 
{{t - u){s - v){to - uo){so - vo)y^ 

\t - to|^"-> - no|2"-2|. - .oP^-> - vo\'^-' (15) 
is finite. This is proved in the following lemma. 

■ 

Remark 1 In the above statement we can replace the normalization j"-2Q+j-e^-2/3+j-e 
6yr-2"+4 5-2/3+4 /(T, 5) where f{T,S) is a deterministic function which converges to zero 
as T, S —)• oo. This is a consequence of the proof below. 

Lemma 1 Let I be given by [75j) . Then for | < a, /3 < | the integral I is finite. 
Proof: Consider the integral 

_i _ 1 

dvQ {{t - u){s - v)) ^ {{to - uo){so - vo)) 4 

X It - toP"-> - noP"-2|s - soP^-> - vo\'^-' 







((t _ _ y)(^to - uo))--^ \t - to\'''-^\n - noP"-2|s - so\^^-' 
X dvo\v - vol"^^'"^ {so - vo)~* . 



Jo 

If sq < V we have 

dvo\v - vo\'^^~'^ {so - vo)~^ = / dvo{v - vof'^~'^ {so- vo)~^ 
Jo 

making the change of variables z = we get 



° f V - so\'^^ l' z{v - So) \ * / V - So 
J \ l-z J \{l-zf 



1-z 

<{v-sof^^^^ 



dz 



I' z-kl - zr^~''dz = {v- sof'-'i /3 ^ - 2/3) 



where /3 is the Beta function. The expresion above is finite for /3 < |. 
Now, if f < So 



so 







dvo\v - wol^^ ^ (so - vo) '4 
dvo {v - vof^'"^ (so - fo)"4 + / dvo {vo - vy^'^ (so - fo) 
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For the first integral in the rieht hand side we make the change of variables z = ^ ^" , and 
for the second one we make z = Then we e:et 

= (so - vf^-l (2/3 -1,^-2/3)+ 2F,{0, 3/4, 2/3 - 1/4; 0)^ 

= (so-z;)2^-t (^/3 (^2/3 -1,^-2/3^ +/3(^^, 2/3-1^ 2^1(0, 3/4, 2/3 - 1/4; 0)) 

wich is finite for ^ < a,/3 < |. Here T is the Gamma function, 2F1 is the Hypergeometric 
function, and we have make use of the property /3(x, y) = ^^(l+^y) formula (see [19], 

formula 1.6.7) 

Vie) f°° 

2Fi(a,6,c;l-x) = ^,,.^^, ^ / w'-\l + wr~'{l + wx)-''dw. 
T{b)T{c-b) Jq 

Proceeding in a similar way with the other seven integrals we conclude that I is 

finite. 



Remark 2 We may note that the prescence of the cos function in the asymptotic behavior 
of the Bessel function does not allow to obtain a renormalization for Ft^s terms of the 
powers of T and S like in JV^. 

The following proposition is a consequence of the Proposition [1] and of the hyper- 
contractivity property of multiple stochastic integrals. 

Proposition 2 For every e > and for ^ < a, /3 < | the sequence 

r-T rS 



1 



^^4a-i+e^4/3-i+e Jq Jq 

converges to zero a.s. as T,S ^ 00. 



Xt,sdBlf 



Proof: As in ^ we can replace the couple {T,S) by a discrete sequence {Tm,Sn) such 
that Tm,Sn- converge to infinity as M,N — )• 00. This is possible since the nominator and 
denominator in ([3]) are continuous a.s. with respect to (T,S). Indeed, the fact that the 
divergence integral in the nominator is continuous follows from [16j, page 293 since the 
integrand X is regular enough, and the integral dtds in the nominator is clearly continuous 
a.s. with respect to the couple (T, 5). For simplicity, we will assume that {Tm,Sn) = 
(M, A^). Let us show that 

Am,n ■■= M-^^+i-iV-^/'+i- / / Xt,sdB^f 

Jo Jo 
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converges to zero a.s. as M, N tend to infinity. We will use the Borel-Cantelli lemma. To 
do this, we will estimate P{Am,n > {MN)~'') for some 7 > 0. For every p > 1 we have 

P{Am,n > (MN)-^) < (MNrBlAMM" 

and since Am,n is a multiple integral in the second Wiener chaos, the inequality ([5]) and 
Proposition [T] implies that 

E|AA/,^r < c{p) {EAlj^^f^ < c{p,a,P){MNy'P. 
Putting together the two above bounds, we get 

^ PiAM,N>{MN)-^)<cip,a,l3) ^ (MiV)P(T-^) 

M>Mo,Af>Afo M>Mo,N>No 
and this series is convergent when 

(e — 7)p > 1 or equivalently 7 < e . 

P 

For every given e > and for p large enough we can always chose a real number 7 such 
that < 7 < e — ^. This, together with the Borel-Cantelli lemma allows us to finish the 
proof. ■ 

The next step is to analyze the denominator in formula ([3]). We have the following 
estimate. 

Proposition 3 For any e > and for any a,f3^ (^, |) 

-, — E / dt I dsX? „ s^oo 00. 

Proof: By the isometry of multiple integrals (jH), the equation ([7]) and the expression of 
the solution to ([I]) 

pt pt ps ps 

EX^^j, = c(a, /?) / du duo dv dvo 
Jo Jo Jo Jo 

xJo{2y/9{t - u){s - v))M2^e{t - uo){s - vo))\u - uol^'^-^lv - ^;o|2^-^ 
making the change of variables i = ^, s = ^ and similar for the other variables we can write 

rT f'S f'l f'l ft ft fs fs 

' / dsdtBXl^ = c(a, /3)r2°+i52/3+i / du duQ dv dvo 

Jo ' Jo Jo Jo Jo Jo Jo 

xJo{2^/9{t - u){s - v)TS)Jo{2^e{t - uo){s - vo)TS) 

x|n-7ioP°->-^o|'^-' 
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and thus 



jr[ dsdtExf, . ^ /"^ /"^ /■* /■* r r 

i = c(a,B) / dt / ds du dun dv dvn 



,h{2^e{t -u){s-v)TS) Jo {2^e{t-u^){s-VQ)TS) 
X 



4 2 



(T5)-4-2 (r5)- 

x|n-noP"->-t;o|'^"'. 
We will use the same idea as in the proof of Proposition [TJ We notice first that 



jQ{2yJe{t-u){s-v)TS) 
(T5)-i-f 

converges to infinity as T, S" — oo by using (fTTI) ; also for T, S large enough 



Jo {2yJe{t-u){s-v)TS) 



iTsy 



is bounded by c{\t — u\\s — v\) ^ almost everywhere s, t, u, w. Also we note that the integral 

1 p\ ft ft fS fS 

dt I ds I du I duo dv dvo\u — uo\'^°'~'^\v — vqI"^^'"^ 
Jo Jo Jo Jo Jo 
_ 1 _i 
{\t-u\\s-v\) 4 (|i - uqIIs - uqI) 4 

is finite for a, /3 S (^j |) by using the same computations as in the proof of Lemma [TJ 
Using Fatou's lemma (we use the following version of the Fatou's lemma: if /„ is a sequence 
of functions such that > —g where g is positive and integrable, then lim J fn > J lim f„) 
this implies that 

J^J^dsdtBXl^ /■! fi /•* /■* r 

lim > c(a, B) dt ds du dun dv dvn 

T,5Too r2"+|-^52/3+|-- - ^ '^'Jo Jo Jo Jo Jo Jo 



j.^ Jo{2^/e{t - u){s - v)TS) Jo{2,Je{t - uq){s - vo)TS) 



1 E .„„^ 1 £ 



T,5^oo {TS)-4-2 {TS)—i-2 



I i2q-2i \2B-2 

X \u — Uq\ \v — Vol = OO. 



At this point we will need the following auxiliary lemma. 

Lemma 2 Consider a sequence of random variables {An)p^ such that P{An > cN^^) < 
OO for some 7 > and for every c > (which implies A]\f ^ almost surely as N ^ 00.) 
Also consider a sequence of a.s. strictly positive random variables {Bi\j)i\i such that ESjv 
00 as N ^ 00. Then 

— yN-^oo almost surely 

Bn 
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Proof: We will use again the Borel-Cantelli lemma. Let C > be arbitrary. Then 



N ^ ' N 

Using the fact that Ei?^ — )• oo as — t- cxd we obtain that 

for N large enoug h. By assumption, Ejv ^ ("^ > N'^,Bn > < oo. Therefore 
and therefore the conclusion follows. 



Remark 3 Lemma \M can be extended without difficulty to two-parameter sequences. That 
is, if Am,n is a sequence of random variables such that YIm N>o Pi^M,N > c{MN)~'^) < oo 
and Bm,n is a sequence of positive random variables such that 'EBm,n -^m,n^oo oo then 

J^^'^ — ?• M N~^oo almost surely. 
Bm,n 

Let us state the main result of this section. 

Theorem 1 Let 9t,s be the estimator given by Suppose that a,/? E |). Then 6t,s 
is a strongly consistent estimator for the parameter 0, that is, 

0T,s ^ almost surely as T, S ^ oo. 

Proof: From relation ([3]), the difference between the estimator and the true parameter is 

^-2a+l^eg-2^+l-e jT jS x.^.dB^^^ 



^ T-'-^\-^S-^^^\-^'}fdt!^dsXl 



a,0 
t,s 



The result is obtained by using Proposition [21 Proposition [3] and Lemma[2] (and the remark 
that follows after this lemma). ■ 
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5 Asymptotic non-normality of the estimator 



We proved in the previous section that the least square estimator 9t,s is strongly consistent. 
This property has been proved in the one-dimensional case in ([9j). Nevertheless, we show 
in this paragraph that the limiting distribution of the estimator is not the same in the one- 
parameter and two-parameter cases. This different behavior is somehow expected since the 
fractional Ornstein-Uhlenbeck sheet does not keep the properties of the fractional Ornstein- 
Uhlenbeck process (for example, the kernel / given by Q can take any real value while 
in the one parameter case it is positive since it is given by an exponential function). In 
order to notice the difference between the one-parameter and the two-parameter case, we 
will focus only on the nominator in the right hand side of ([3]). To check the asymptotic 
normality we use the following criterium: (see Theorem 4 in |17j). 

Theorem 2 Let {F^^k > 1), Fj. = Inifk) (with G 7^®" for every k > 1) be a sequence 
of square integrable random variables in the n th Wiener chaos such that E[F^] 1 as 
/c — >• oo. Then the following are equivalent: 

i) The sequence (-Ffe)A;>o converges in distribution to the normal law M (0,1). 
a) converges to n in Lp'iyt) as k ^ oo. 

Denote by 

Proposition 4 Assume ^ < a, /3 < | and let Ft^s be given by [W\). Then when T, S tend 
to infinity, 

1 

Ft,s 

does not converges in distribution to the normal law A^(0, 1). 

Proof: We will prove that \\D -^^FT^s\\\a,fi does not converges to 2 in L^(r2) as T, 5 — )• oo. 
Let us denote by / the symmetrization of / with respect to the variables (t, s), (u, v) 

f{{t,s), {u,v)) = ^ {f{t,s,u,v) + f{u,v,t,s)) . 

We have 

Dt,s^FT,s = ^Dt,sh{f) = ^Dt,sh(f) =^{h{f{;;t,s))+h{t,s,;-)). 
Here "•, •" represents the variable with respect to which the integral Ii acts. 
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We obtain 



1 1 1'^ 1'^ 

\D FT,s\\\ic,f} = — 9 — c(a)c(/?) / dt I ds dto / dso 

^T,s cft.s Jo Jo Jo Jo 



'T,S 

X [h (/(•, •, t, s)) + h {fit, S, ; •))] [h ifi; ■,to, so)) + h if {to, SO, ; ■))] 

x|t-to|2«-2|^_5o|2/3-2 



r rS f-T fS 



c(a)c(/3) / dt I ds I dto I dso 



^T,s Jo Jo Jo Jo 

X [I2 (/(•, •, t, s) /(•, ■,to, so)) + h ifi-, ■,t, s) /(to, so, •, •)) 

+h if it, S, ; •) fi; ■,to, so)) + h if it, S, ; ■) /(to, SO, ■))] 

x\t- top"-2|s - so\^^-^ + E\\D^Ft,s\\hc,0 

Let us note that that ^WD-^^FT^sWy^a,!} is equal to 2. This follows from the fact that for 
any multiple integral of order n we have 

B\\DlM)\\n = nBInif)^. 

It remains to show that the terms containing multiple integrals of order 2 does not 
converges to zero in L'^iQ) as T, S" — t- oo. The first and the fourth summand are similar, as 
they are the second one and the third one. We will handle only the first summand denoted 
by ^2^5 (because the other three terms can be studied analogously). We can write, using 

the definition of the scalar product in the Hilbert space (^"''5)®^ and the expression ([10]) 
of the kernel / in terms of the Bessel function Jo 
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[ dt [ ds [ dto [ dso\t-tof''-^\s- sof^-"^ 
Jo Jo Jo Jo 



'T,S 

l2ifi;-,t,s)(S>fi;-,to,so)) f 



rT rS rT rS rT pS rT pS 

2 dt ds dto / dso / du dv I duo / dvo 
Jo Jo Jo Jo Jo Jo Jo Jo 

(/(•, ■,t,s) (g) f{-,-,to,So),f{-,-,U,v) f(-,-,Uo,Vo))^^^ g^02 

rT rS rT rS rT /vS' rT rS 



C— — / dt ds dto / dsQ / du dv duo I dvo 
'^T,s Jo Jo Jo Jo Jo Jo Jo Jo 

pT pS pT pS pT pS pT pS 

/ dx dy dxo / dyo / da db dao / dbo 
Jo Jo Jo Jo Jo Jo Jo Jo 

f{x, y, t, s)f{xo, yo, to, so)f{a, b, u, v)f{ao, bo, to, so) 

\t - to\'--'\s - so\'^-'\n - ^o|^"-> - vo\'^-' 



\x — a 



1 



,2a-2 



\y-b\^^-'\xo-ao\"^-'\yo-bo\'^~' 



T rS rT rS rT rS rT rS 



- C— — / dx dy dxo / dyo / da db dao / dbo 
^T,S Jo Jo Jo Jo Jo Jo Jo Jo 

i-x i-y i-XQ i-yo ra rb r-aa rbo 

/ dt ds dto / dso / du / dv / d.uo / dvo 
Jo Jo Jo Jo Jo Jo Jo Jo 

Jo{2 ^e{^ t){y - s)) Jo{2^ e{xo - to){yo - so)) 

Jo (2 ^d{a-u){b-v)) Jo {2^e{ao - uo){bo - vo)) 

\t - tor-'\s - so\'^-'\u - nop"-> - vol'f"-' 

= ^ and similar for the other variables we get 

J_j.8a^8/3 f f f f f f j f 

^T,S Jo Jo Jo Jo Jo Jo Jo Jo 

rx ry nxo ryo ra rb /-ao rbo 

/ dt ds dto / dso / du dv duo / dvo 
Jo Jo Jo Jo Jo Jo Jo Jo 

Jo {2y /9{x - t){y - s)TS) Jo {2^ ejxQ - to){yo - so)TS) 
Jo{2VO{a - u){b - v)TS)Jo{2y/0{ao - uo){bo - vo)TS) 
\t - to\"^-'\s - so\'^-'\u - uo\'-'\v - vo\'^-' 
\x - a|2"-2|y - 6|2/3-Vo - aor-'\yo - bo\'^-' 
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Using the asymptotic behavior of the Bessel function when its variable is close to infinity 
(see (fTTI) ) we see that 



E 



A. 



(1) 

T,S 



S Jo Jo Jo 



dxo / dyo / da db dao / dbo 
Jo Jo Jo Jo Jo 

ex i-y j-X(i f-yo r-a r-b i-ao i-bo 

dt ds dto / dso / du dv duo / dv^ 
'0 Jo Jo Jo Jo Jo Jo Jo 

-1 

((x -t){y - s){xo - to)(yo - so){a - u){b - v){ao - no)(6o - ^o)) * 

U J. \2a-2\ \20-2\ \2a-2\ \2B-2 



\x — a 



\2a-2 



i.\2S-2\ \2a-2\ i \2B-2 

y-b\P ko-aol |yo - ^o| ^ ■ 



and considering the fact that the last integral is finite for ^ < a, /3 < | (the proof of this 
fact is similar to the proof of Lemma [T] ) it is straightforward to see that, for T, S close to 



infinity, the quantity E 



A. 



(1) 

T,S 



does not converges to zero. 



Acknowledgement: The authors would like to thank Prof. Soledad Torres for 
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